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In this article, we present a method for detecting the amount of polymer added in the biomass pyrolysis 
process to obtain fuel for energy generation. The main aim of these research was to develop a new 
methodology to determine the amount of polypropylene added in the pyrolysis process, which was used 
as a case study of the type of impurities that can contaminate biomass. In our best knowledge these 
methods, based on chemometric methods combined with spectral data from Fourier transform infrared 
spectroscopy (FTIR) applied to bio-char analysis are new one and could allow a more accurate control of 
bio-char than existing methods such as selective dissolution analysis of the 14 C content. As a chemomteric 
methods in this studies the classification analysis (CA) and principal component analysis (PCA) were 
applied. Thereby we could describe relation between the third principal component (calculated from 
FTIR spectra) and amount of polypropylene additive depending on the pyrolysis temperature. 

To the best of our knowledge, there have been no reports in the literature on determining the amount 
of polymer added in using a bio-char control instead of raw biomass for co-firing. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

The global primary energy demand in 2013 increased by 2.3% 
over the previous year. Market analysts predict that the global 
energy demand will increase by 41% by 2035. Rising fossil fuel 
prices have led to the expectation that by 2035, renewable energy 
sources will become the second largest source of electricity 
generation globally [1-3]. The increasing use of biomass for the 
production of primary energy has resulted in more frequent 
problems that are related to the physicochemical properties of 
biomass. These problems may be eliminated (or at least limited) 
by pre-processing the biomass into fuels with greater stability, 
as characterized by enhanced physicochemical properties (e.g., 
the calorific value). There are well known methods for thermal 
biomass upgrading (such as torrefaction or pyrolysis). These pro¬ 
cesses can be used to obtain fuels with quality parameters that are 
comparable to those of hard coal. These parameters depend on the 
process conditions and the properties of the material used. The key 
to the wide application of bio-chars as alternative renewable fuels 
for energy generation is the type of material used in thermal pre¬ 
processing. Bio-char must be completely produced from biomass 
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to be considered as a renewable fuel. Unfortunately, the bio-char or 
torrefied biomass entering the market may be knowingly contami¬ 
nated by chars obtained from coals or plastics pyrolysis. Therefore, 
an appropriate (and relatively cheap) method of analysis must be 
developed to accurately classify the bio-material used for energy 
production. 

The first mention about biomass origin analysis in which the 
statistical and chemometric methods were applied have been pre¬ 
sented in the work on the biomass properties [4,5 . Unfortunately 
authors used chemometric methods only to data visualization, do 
not builded a specific algorithm for classification, which would 
determine the type of biomass on the basis of chemical analysis 
of samples. In other manuscripts can be found some information 
about spectrometry method coupled with chemometric analysis 
which were used to quality analysis of biomass [6-9]. Nowadays 
the analytical laboratories that perform analysis of biomass and 
biodegradable materials derived from biomass have the ability to 
mark the biodegradability, which is determined according to the 
European standard: EN 15440 Testing for Solid Recovered Fuels 
(Selective Dissolution Method (SDM) and Carbon-14 Method). 
According to the standard, this method is not appropriate if the SRF 
sample has biomass components that are insoluble in sulfuric acid 
or fossil-based components that are soluble in the acid. The sec¬ 
ond method - the Carbon-14 Method - measures the radiocarbon 
content of the mixed wastes and is applicable to all materials [10]. 
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Table 1 

Experimental matrix (central composite rotatable design 2 2 + 2 x 2 + 3) for biomass 
and biomass-PP blends. 



No. 

Coding - classical 
experimental design 

Temp. [°C] 

PP [%] 


1 

500 

10 


2 

500 

50 


3 

700 

10 


4 

700 

50 

Regular experiment 

9 

458 

30 


6 

740 

30 


7 

600 

0 


8 

600 

58 


9 

600 

30 

Central 

10 

600 

30 


11 

600 

30 


12 

500 

0 


13 

600 

0 

Blank 

14 

700 

0 


15 

458 

0 


The analysis methods that were initially used were based on 
chemometric analysis, and the pre-classification of the biomass 
and the origin of the bio-char were performed using proximate and 
ultimate analysis and spectroscopy [11-14]. 

In this paper, we present current results from using chemomet¬ 
ric analysis to classify and verify the origin of bio-char: a FT-IR DRIFT 
analysis was conducted in combination with a few chemometric 
methods. 

2. Experimental 

2.1. Materials and analytical methods 

In this study on the pyrolysis process, the materials used 
consisted of pine wood biomass and polypropylene without any 
filler, in the form of 4 mm radius granules and colored with char¬ 
coal (Daplen provided by Borealis). The biomass investigated was 
prepared for chemical analysis and the thermal conversion pro¬ 
cess. A detailed description of the methods used in the raw material 
preparation has been previously presented [15]. 

The prepared samples were pyrolyzed under different ther¬ 
mal conditions ranging from 458 °C to 740 °C and with different 
amounts of added polypropylene, following a central composite 
design in the design of experiment method (DOE, see Table 1). In 
the thermal conversion process (co-pyrolysis), 150-g samples were 
pyrolyzed under a nitrogen atmosphere. A detailed description of 
the pyrolysis process, a schematic of the laboratory-scale apparatus 
and the thermogravimetric analysis of the applied polypropylene 
has been given in a previous paper [16]. Table 2 summarize the 
primary parameters used in the design and the subsequent perfor¬ 
mance of the experiments. 

A similar research study on biomass pyrolysis under various 
thermal conditions has already been published [17-19]. Flowever, 
only pure biomass without the addition of any other materials was 
investigated in the aforementioned study. All of the samples of 
carbonized biomass were investigated by FT-IR analysis using the 
DRIFT method. 

2.2. Analytical method for bio-char investigation 

A sample (0.6 mg) was mixed with KBr (300 mg), finely ground 
and placed in a sample cup. The sample cup was overfilled with 
the sample, and the sample surface was flattened with a spatula. 


Table 2 

Principal component analysis results. Significant values of the principal components 
(above 0.7) are marked in bold. 


No 

PCI 

PC 2 

PC 3 

Temp. [°C] 

PP I%] 

1 

0.984 

0.141 

-0.079 

500.00 

10.00 

2 

0.920 

0.374 

0.086 

500.00 

50.00 

3 

0.933 

0.262 

-0.035 

458.00 

30.00 

4 

0.978 

0.164 

-0.106 

458.00 

0.00 

5 

0.989 

0.090 

-0.114 

500.00 

0.00 

6 

0.613 

0.746 

0.183 

700.00 

0.00 

7 

0.453 

0.853 

0.242 

600.00 

30.00 

8 

- 0.111 

0.744 

0.588 

700.00 

10.00 

9 

-0.020 

0.964 

0.225 

600.00 

58.00 

10 

0.501 

0.836 

0.199 

600.00 

30.00 

11 

0.496 

0.850 

0.141 

600.00 

0.00 

12 

-0.050 

0.233 

0.962 

740.00 

30.00 

13 

-0.078 

0.264 

0.959 

740.00 

30.00 


The sample was then placed into a diffuse reflectance accessory 
to generate the spectrum using potassium bromide as a refer¬ 
ence material. The relationship between the concentration and the 
spectral response was linearized by applying the Kubelka-Munk 
function to the spectral data. Fig. 1 presents a representative 
spectrogram that was used to analyze the carbonized biomass con¬ 
taining PP (500 °C - 0% PP, 10% PP, and 50% PP). 

2.3. Chemometric and compositional data analysis of bio-char 

First, all of the FT-IR spectral data were normalized using Eq. (1 ), 
which was essential for the subsequent chemometric analysis and 
to properly interpret the results [20]. This procedure was used to 
unify examined variable values and determine the variance in the 
data set considered. 


The terms in the equation above are defined below: 


• Xy - standardized parameter value, 

• a 2 j- - initial value of parameter, 

• bj - average value of parameter, and 

• Sj - standard deviation of j-th parameter. 


Fig. 2 shows the results of the spectral data normalization that 
was performed using Eq. (1). Fig. 2a shows the raw spectral data 
(before normalization) for the tested samples, and Fig. 2b shows the 
spectral data that were obtained after the normalization process 
and were then subjected to cluster analysis and principal com¬ 
ponent analysis. The normalization procedure enabled the tested 
materials to be analyzed comparatively using the raw spectra from 
the FT-IR analysis. 

After pre-preparation, the data set was subjected to cluster anal¬ 
ysis (CA) followed by principal component analysis (PCA). 

Cluster analysis is one of several exploratory multivariate 
methods that can be used to determine the relationships among 
variables. A data set must be normalized to examine the similarity 
between variables. In hierarchical clustering analysis (HCA), Ward’s 
method is used to obtain a cleaner plot of clusters. This method dif¬ 
fers from other methods in that an analysis-of-variance is used to 
evaluate the distance between clusters. In this method, the sum of 
squares (SS) of any two clusters that can be formed at each step 
is minimized. This method is considered to be very efficient. The 
Euclidean distance given by Eq. (2) was used in the analysis: 


dij — 


\ 




k =1 


( 2 ) 
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FT-IR DRIFT 



Wavenumber [cm' 1 ] 


— 500/10 

— 500/50 

— 500/0 


Fig. 1. Representative spectrogram of char from biomass and polypropylene. 
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where 

• m - number of variables, 

• dij - Euclidean distance between objects i and j, 

• - value of variable i, and 

• Xj k - value of variable j. 

It is quite difficult to directly analyze a large data set (e.g., a 
spectral data set) is quite for a large number of test materials. Hence, 
the interpretation of the data that are obtained from cluster analysis 
was simplified by presenting the data in the form of dendrograms, 
as shown in Fig. 3, instead of in tables. 

The smallest value is the Euclidean distance that corresponds to 
the largest distance between two objects. The Euclidean distance 
can be graphically interpreted as the length of the vector starting 
at i and ending at j. 

The second technique that was applied in this study was princi¬ 
pal component analysis, which was found to be a very efficient tool 


for multidimensional data analysis. In brief, PCA is a data reduc¬ 
tion technique in which new variables are calculated from linear 
combinations of the original variables. In the PCA technique, the 
eigenvalue and eigenvectors are extracted from the covariance 
matrix of the original variables. The salient characteristic of PCA 
is its ability to determine the primary components (i.e., the specific 
values) while decreasing the information carried by subsequent 
variables. In some cases, a new variable may take on a specific phys¬ 
ical or chemical interpretation as a result of PCA analysis [21,22]. 

Principal component analysis can be used for 


• the dimensional reduction of data space (i.e., to reduce the num¬ 
ber of variables), 

• the creation of new transformed variables, (i.e., the principal com¬ 
ponents), and 

• the graphical presentation of multidimensional data in a 3D sur¬ 
face without significant loss of information. 


FT-IR DRIFT after normalization 
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Fig. 2. Spectra of investigated materials (bio-char) (a) before and (b) after normalization. 
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Fig. 3. Representative plot of horizontal hierarchical tree. 


3. Results and discussion 

The data from the FT-IR analysis of the investigated carbonized 
biomass were normalized after being subjected to cluster analysis 
followed by principal component analysis. All of the chemometric 
and statistic operations described in this paper were carried out 
using STATISTICA software [23, 24]. 

After the cluster analysis, three data subsets emerged that 
were characterized by different temperatures in the process. The 
results were presented in the form of a tree diagram: Fig. 4 shows 
a dendrogram instead of the common tabular form. The figure 
shows that the first group (shown in blue) contained samples that 
were obtained during pyrolysis at temperatures up to 500 °C. The 
second group shown in the diagram (in green) describes the sam¬ 
ples obtained from pyrolysis at a temperature of approximately 
600 °C. Significant differences in the distances between the tested 
samples originated from the use of different PP contents in the 


biomass mixture. The unstable conditions during pyrolysis cannot 
be neglected (shown in red in Fig. 4). The final and third group 
(shown in red in Fig. 4) consisted of the samples that were obtained 
during pyrolysis at temperatures above 700 °C, with one exception. 
This exception was a biomass sample with a 58% PP content that 
was pyrolyzed at 600 °C. This high polymer content changed the 
process conditions, resulting in a larger amount of processed 
material. FICA was used to make approximate predictions of the 
temperature and the added amount of PP for the new samples that 
obtained in pyrolysis process and subjected to FT-IR analysis. 

The polymer that was used in this part of the study (Borealis’ 
Daplen) undergoes depolymerization (thermal decomposition) at 
approximately 400 °C. This process releases free radicals that 
quickly react with cellulose, hemicelluloses or lignin particles that 
are present in the biomass. Dehydration and decarbonization reac¬ 
tions result in the polypropylene radicals lowering the oxygen 
content in the carbonized material. The relationships among the 
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Fig. 4. Horizontal hierarchical tree plot for different types of samples. (For interpretation of the references to color in the citation of figure, the reader is referred to the web 
version of the article.) 
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properties of the bio-char resulting from biomass thermal conver¬ 
sion and the amount of polypropylene added in this process have 
been briefly discussed in other publications [16,24]. 

The spectrum of biomass carbonized together with PP exhib¬ 
ited a similar pattern as the pure biomass that was processed at 
a higher temperature. This observation was confirmed by elemen¬ 
tary analysis (C, H, N, S, O), which demonstrated a clear relationship 
between the PP additive and the oxygen content in the tested mate¬ 
rial. The elemental analysis was conducted using a direct method 
with a Macro Cube CHNS, O, Cl elemental analyzer that was pro¬ 
vided by Elementar Analsensysteme GmgH (for high temperature 
pyrolysis), and the ND-IR method was used for oxygen detection 
[25]. 

Following cluster analysis, the second type of chemometric anal¬ 
ysis, principal component analysis, was performed on the spectral 
data set. This analysis produced a new data set with fewer variables: 
the basic tree variables obtained were the primary components. 



Observed temperature [°C] 

Fig. 7. Scatter plot showing observed temperature values and predicted tempera¬ 
ture values obtained using Eq. (3). 


The data set that was obtained by the analyzing the investigated 
materials using the primary components are presented in Table 2. 

There was a total of 3 PCs in this case, which described 97.89% of 
the total variance of the data set. The principle component analysis 
results confirmed the cluster analysis results. Fig. 5 shows that the 
agglomerates resulting from principal component analysis corre¬ 
sponded to those obtained from cluster analysis and are shown in 
Fig. 4. 

The principle component analysis and the wavenumber from 
the FT-IR analysis showed that the first component contained 
information at 1593 cm -1 (C=C groups), the second principle com¬ 
ponent contained information at 3630 cm -1 (free —OH groups), and 
the third principle component contained information at 3439 cm -1 
(—OH groups). Additionally, Table 2 (and the diagram in Fig. 5) show 
that the new principal component corresponded to the pyrolysis 
conditions under which the samples were obtained. The first prin¬ 
cipal component (PC 1) corresponded to char that was obtained at 
500 °C, PC 2 corresponded to char that was obtained between 600 
and 700 °C, and PC 3 corresponded to char that was obtained below 
700 °C. 

Unfortunately, the obtained values of the PC’s could not be used 
to describe the amount of PP that was added in the pyrolysis pro¬ 
cess: this information could only be obtained using subgroups for 
specific temperature ranges. This result was confirmed using a 


"PC1"*"PC2" 

"PCI" 

"PC2" 

"PC2" A 2 

"PC1" A 2 

"PC2"*"PC3" 

"PCs- 

"PC3" A 2 

"PC1"*"PC3" 

p=,05 

t statistic 

Fig. 6. Pareto charts of (a) significant temperature effects and (b) PP added in pyrolysis process. 
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a) 


Chart for 458-500°C 



b ) Chart for 600°C 



4 Chart for 700-740°C 



Fig. 8. Change in the value of the third major component depending on the amount of PP added in different temperature (a-c) of pyrolysis. 


general regression model (surface response regression) and was 
presented in the form of Pareto charts (Fig. 6). 

Fig. 6a shows that the first two components and their combina¬ 
tions had a significant effect on the predicted temperature at which 
char (or bio-char) was obtained (p > 0.05). Fig. 6b shows that none 
of the designated new variables could be used to determine the 
total amount of added polypropylenes in the pyrolysis process. The 
applied multivariate regression method was used to develop Eq. 
(3) (see below) for the temperature at which the pyrolysis process 


was conducted using only spectral data from the FT-IR analysis and 
the PCA analysis, where R = 0.986 and R 2 = 0.972 (Fig. 7). 

Tpy W = 7304.2 - 10145.6 • PCI + 3468.8 • PCI 2 - 10387.1 • PC2 
+ 3607.8 • PC2 2 - 4316.8 • PC3 - 1518.9 • PC3 2 
+ 8546.9 PCI PC2 +3412.6 PCI • PC3 
+ 4267.9 PC2 PC3 (3) 
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As previously mentioned, the amount of polypropylene that was 
added in pyrolysis could not be directly calculated using the data 
obtained from principal component analysis. When the pyrolysis 
temperature was known (or the value could be estimated using Eq. 
(3)), the third principle component of the three temperature ranges 
(458-500 °C, 600 °C, 700-740 °C) could be used to estimate the 
amount of added PP in the biomass pyrolysis, as shown in Fig. 8a-c. 

4. Conclusion 

In this paper, we have presented preliminary research results 
obtained using chemometric methods to determine the presence 
and amount of polymeric materials that were carbonized together 
with biomass in pyrolytic processes. Cluster analysis was used to 
determine that three well-defined groups of biomass were car¬ 
bonized up to 500 °C, 600 °C and 750 °C. The predicted amount 
of polypropylene that was added to the biomass in the thermal 
conversion process was determined using FT-IR spectral analysis 
with standard normal variate (SNV) transformation as the data pre¬ 
preparation method, coupled with chemometrics methods (HCA 
and PCA) that could not be properly carried out over a wide tem¬ 
perature range. The study results showed that the polypropylene 
amount in char could only be predicted over narrow temperature 
ranges, for example 400-500 °C, 500-600 °C, and 600-700 °C. 

In future studies, we will conduct additional tests on new sam¬ 
ples of chars obtained from the thermal conversion of biomass and 
a few types of plastic waste, such as PET, ABS, and SBR. 
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